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Abstract. Temporal behaviour of GRB afterglow light curve is derived 
for the case where the electron energy distribution is relatively hard, with the 
power-law index p lying between 1.0 and 2.0. It is shown that the expected 
behaviour will be the same as that for p > 2.0 if the upper cutoff in the electron 
energy distribution evolves in direct proportion to the bulk Lorentz factor of 
the blast wave. 
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1. Introduction 

Detailed observations of afterglows of Gamma Ray Bursts over the last four years have 
established that they exhibit power-law broadband spectra and power-law temporal decay 
of their light curve. The generally accepted model for the afterglow, called the fireball 
model, explains this emission as being due to synchrotron emission from a relativistically 
expanding blast wave which accelerates electrons to large Lorentz factors, with a power- 
law energy distribution (see Piran 1999 for a review). Recently in several cases the light 
curve of the afterglow has been seen to undergo a break into a steeper power law, a 
behaviour that is expected if the burst is beamed into a narrow solid angle (see Rhoads 
1999, 2001). Theoretical predictions for the spectral and temporal evolution have been 
made in detail for both isotropic (Wijers, Rees & Meszaros 1997; Waxman 1997a, b,c; 
Sari, Piran & Narayan 1998; Wijers & Galama 1999) and beamed (Rhoads 1997,1999; 
Sari, Piran & Halpern 1999) fireballs, and these have enjoyed considerable success in 
modelling the behaviour of observed afterglows. 
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Nearly all theoretical work in the literature so far assume that the energy distribution 
of the injected electrons is a power-law: 

^(7e) OC -f-P, (7m < 7e < 7«) (1) 

(where 7e is the Lorentz factor of the electron and 7„j and 7„ are the lower and upper 
cutoff of the energy distribution respectively), with an index p larger than 2.0. This 
assumption simplifies the derivations, since the particles at the low-energy end domi- 
nate both the total number and the energy content for such a steep energy distribution. 
The values of p derived from observations of most afterglows do indeed fall above 2.0, 
thereby allowing meaningful comparison being made between theoretical predictions and 
observations in these cases. 

However, at present no compelling argument is known as to why the energy distribu- 
tion of the accelerated electrons must always be so steep. Indeed a fairly large dispersion 
is seen in the spectral index distribution of shock-accelerated electrons in Galactic shell 
supernova remnants, which include several cases where p is inferred to be less than 2.0 
(cf. Green 2000). Moreover, in Crab-like nebulae, where the particle energy distribution 
is shaped possibly by a relativistic standing shock (Rees & Gunn 1974, Kennel & Coro- 
niti 1984), the value of p is almost always found to be less than 2.0. In the context of 
GRB afterglows, p ~ 1.5 has been invoked for GRB 000301c (Panaitescu 2001) and for 
GRB 010222 (Sagar et al 2001, Cowsik et al 2001). 

Clearly, modelling of such a hard spectrum afterglow at present suffers from the 
handicap that theoretical predictions specific to such energy distributions are not available 
in the literature. Panaitesu (2001) makes a detailed case study of GRB 000301c with 
p ~ 1.5, but does not provide general results that are easily applicable to other cases. The 
aim of this paper is therefore to extend the predictions of the fireball model to the case of 
p < 2.0. In this paper I will address only the most commonly used spectral and dynamical 
regimes, namely slow cooling, adiabatic evolution for both isotropic and beamed bursts. 
I will consider a range of p between 1.0 and 2.0 and make certain simplifying assumptions 
that allow easy analytical treatment. I will also assume that the afterglow is optically thin 
over the entire range of frequencies of interest. A more detailed and exhaustive study of 
hard spectrum fireball models will be reported in a future publication (D. Bhattacharya & 
K.M. Basu, in preparation). In what follows I will adopt, wherever applicable, expressions 
from Rhoads (1999) and Wijers & Galama (1999) with slight modification in notation. 



2. Adopted results 

I list below the expressions already available in the literature which I adopt for the 
purposes of the present paper. These expressions are not affected by the change of 
energy distribution index p. 
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2.1 Dynamics 



I will confine myself to the adiabatic regime of expansion of the blast wave. Rhoads (1999) 
presents results for dynamical evolution both before and after the light curve break, which 
corresponds to a time tb (in the frame of the observer) when the initially tightly collimated 
ejecta begins to expand predominantly sideways. The expressions corresponding to times 
before tb can also be used to represent isotropic bursts by setting the initial solid angle 
of coUimation to 47r. 



From Rhoads (1999) I adopt the following expressions for the dynamical evolution of 
the blast wave: 



r = 2-5/4 
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where T is the bulk Lorentz factor of the blast wave, Tb the value of T kit t = tb, Eq is the 
total energy of the blast wave, is the initial opening angle of the collimated ejecta, p 
is the ambient density and Cg is the sound speed of the postshock medium, also taken to 
be the speed of lateral expansion, c, as usual, is the speed of light. For a spherical blast 
wave, the appropriate expressions can be obtained by substituting ttOq with An. The time 
t is measured in the frame of the earthbound observer, z is the redshift of the afterglow. 



2.2 Magnetic Field 



Using the expressions in Rhoads (1999) the evolution of the postshock magnetic field in 
the blast wave (as measured in the comoving frame) can be written as: 

/ c \ 1/2 
B = f STT^espj Tc (6) 

where es is the fraction of the postshock thermal energy converted into magnetic energy. 
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2.3 Radiation 



As explained by Wijers & Galama (1999), the observed location of the peak of the syn- 
chrotron spectrum radiated by a single electron of Lorentz factor 7e is 

0.286 e 2 

^{le) = — rB7^ (7) 

1 + Z TTTTleC 

Integrated over the power-law energy distribution of electrons, one obtains a power-law 
radiation spectrum for the whole afterglow, with the peak lying at 

ri?7f„ (8) 



1 + Z TTirieC 

where jm is the lower cutoff of the energy distribution (see eq. (1)). The factor Xp is a 
function of the index p of the energy distribution. For 1.0 < p < 2.0 the value of Xp lies 
between ~ 2.0 and ~ 0.65 (Wijers & Galama 1999). Here e and rrie are the charge and 
the mass of the electron respectively. The received flux per unit frequency at this peak 
of the afterglow spectrum is given by 

where Ng is the total number of radiating electrons, 11 is the solid angle in which the 
radiation is beamed and d is the luminosity distance to the afterglow from the observer. 
(j)p is a p-dependent factor, and lies between 0.4 and 0.6 for 1.0 < p < 2.0 (Wijers & 
Galama 1999). According to Rhoads (1999) Fm works out to be, for t < t^, 

where fie is the mean molecular weight of the ambient medium and rup is the proton 
mass. Additional factors of order unity will need to be inserted in eqs. (§)-(l^ to 
represent the result of integration over different parts of the fireball. Eq. ( p^ shows that 
Fm is independent of time, i.e. the flux at the peak is constant. If we denote this value 
as F^, then the evolution after tb can be written as (Rhoads 1999) 

(11) 



The electron energy above which synchrotron cooling is important within the expan- 
sion time corresponds to the Lorentz factor (Wijers & Galama 1999) 



where ctt is the Thomson scattering cross section. 
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Using eq. (|^) and the expressions of the comoving magnetic field given above one 
obtains the expression for the coohng frequency from 



0.286 X 384ci/2 /c^\3/2 



2?ne /Cs Y 

(l + z)5/2(40)3/2 4 \~) 

[t < h) (13) 
>^c{tb) — constant {t > tij). (14) 



-3/2 t^O 1/2 



3. Results for Flat Spectral Index 

We now have all the pieces necessary to compute the evolution of the afterglow spectrum 
for a hard (p < 2) energy distribution of electrons. For 1 < p < 2 and 7„ 3> 7m 
the only way this modifies the evolution is by changing the evolution of 7m with time. 
As in Sari, Piran & Narayan (1998) we note that the postshock particle density and 
energy density are 4Tn and AT^nmpC^ respectively, where n is the number density of 
the ambient medium. Assuming a fraction eg of the postshock thermal energy goes into 
power-law electrons, these quantities can be equated to integrals over the electron energy 
distribution: 

iV(7e)d7e = 4rn (15) 

7u 

7eTOeC^Af(7e)d7e = ee'^^^nrUpC^ (16) 

7m 

Clearly, for 1 < p < 2 the dominating limit in the first integral is 7m while that in the 
second integral is 7„. Using the fact that 7m ^ 7u one then obtains 



7r> 



p-lj me 



(17) 



This is the key element that causes the difference of evolution between the hard spectrum 
and steep spectrum afterglows. To recall (e.g. from Sari, Piran & Narayan 1998), for 
p > 2, the value of 7m evolves as 

7™=.e(^)^r (18) 

The integral in eq. ( [l6| ) is carried out over the injected energy spectrum of electrons, 
which is an unbroken power-law up to 711. This, therefore, is a measure of the total 
energy the acceleration process injects into relativistic electrons, and for the purposes of 
this paper I assume that this is a constant fraction (ee) of the postshock thermal energy. 
The spectrum of the accumulated electrons, however, would steepen to 7e ^^"^^^ beyond 
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the coohng break 7c (eq. |T^) because of the radiation losses suffered after acceleration. 
Over the period of interest, 7c would in general be much less than 7„ (see, e.g. Gallant 
& Acherberg 1999, Gallant, Achterberg & Kirk 1999), so for p < 2, only a small fraction 
of the total injected energy will remain in the accumulated electrons. Dai and Cheng 
(2001) have computed the evolution of the afterglow assuming that the ratio (ec) of this 
remaining energy to the postshock thermal energy stays constant with time. While the 
constancy of either Ce or Ec as defined above is a questionable assumption, the degree of 
difhculty in arranging a physical situation to maintain a constant is certainly greater. 
We therefore derive our results assuming a constant ee, although the final results will be 
general enough for application to either case. 

Eq. ( [l7| ) shows that the evolution of 7m in the hard spectrum case depends on how 
7„ changes with time. This depends on the details of the particle acceleration process in 
the ultrarelativistic blast wave, which have so far not been very well understood (see the 
review by Bhattacharjee and Sigl (2000) and references therein). Broadly speaking, the 
maximum energy achieved by an electron in the acceleration process would be limited 
either by radiation losses within the acceleration cycle time or by the cycle time itself 
exceeding the age of the blast wave. These quantities depend on the shock parameters 
as well as the upstream magnetic field strength. Since the evolution of most of the shock 
parameters can be expressed as power-law dependences on F, for the purposes of this 
paper I make the simplifying, but perhaps not very unreasonable assumption that 7,j for 
a given afterglow is a function of T alone, and parametrize this dependence as a power 
law: 

In = er« (19) 

where ^ is a constant of proportionality. The value of q, however, may not be constant with 
time, and may depend on the dynamical regime. For example, in the simplest acceleration 
models (cf. Gallant and Achterberg 1999), q ^ 0.5, independent of dynamical regime, if 
the acceleration is limited by radiative losses; but if the age t of the blast wave limits the 
acceleration, 7„ oc Ff, which yields a dynamics dependent q. 



Eq. (|T9|) yields 



p(l+pg-2g)/(j9-l) 



(20) 



The dependence of 7m on F reduces to that for p > 2 (eq. |l^) if g = 1, i.e. if the upper 
cutoff energy also is directly proportional to the bulk Lorentz factor of the shock. In this 
case all the results derived for the temporal behaviour of p > 2 afterglows will also be 
applicable to those with p < 2. 

It is now straightforward to obtain the dependence of Vm on time by inserting eq. ( p0| ) 
in eq. (@), and using the appropriate expressions for F and B. The result is 



1 + Z TT rrie 
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for t <tb and 



1 + Z TT rUeC 



p-1 
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2/(p-l) 



^^p(p+l+2M-4g)/(p-l) 



-(p+P9-2g)/(p-l) 



(22) 



for t > tb- Here stands for B{tb). 

Noticing now that befow and above the coohng break the aftergfow flux is given by 

-(p-l)/2 



F = F 



F = F 



-(p-l)/2 / , X -p/2 



(i/c < 1/ < Uy) 



(23) 

(24) 



(where Vu is v {ju ) ) , wo can easily obtain the time dependence of afterglow flux by inserting 
the time dependence of Fm, i^m and Vc- The final results are 



Oc j,-(p-l)/2i-3(p+p9-2,)/8 (^^<j,<^J 



OC V 



p/2^-[3p+2+3g(p-2)]/8 



for t <tb and 



F^ OC ,,-(P-l)/2f-b+2+g(p-2)]/2 (j.^<z.<j.^) 



(25) 
(26) 



(27) 
(28) 



for t > tb- As one can verify, these expressions reduce to the familiar expressions for p > 2 

by setting q = 1. Further, in the case of constant Ec (Dai and Cheng 2001) 7c plays the 
role of 7i,, and the corresponding results can be obtained by inserting the dependence of 
7c on r in the above equations: q = —1/3 for t < tb and q = —1 for t > tb- 



4. Conclusions 



I have presented above the expected behaviour of GRB afterglow light curves when the 
index p of the power-law energy distribution of electrons lies in the range 1.0 < p < 2.0. 
The results presented here correspond to the optically thin, adiabatic, slow cooling regime. 
The total energy content in such an energy distribution is dominated by the upper cutoff 
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Lorentz factor 7„, and hence the evolution of 7„ influences the evolution of the light 
curve. I derive the light curve behaviour assuming a simple power-law dependence of 7„ 
on the bulk Lorentz factor F of the blast wave. It follows that the behaviour of the light 
curve for 1 < p <2 will be similar to that for p > 2 if 7„ cx T. 

It ought to be remembered that the broken power-law description of the afterglow 
spectrum and light curve presented here represents only the asymptotic behaviour, in 
reality the transitions between different regimes are expected to be smooth. Moreover, 
for relatively hard electron energy distributions considered here, synchrotron cooling is 
expected to cause a pile-up of particles at the cooling break 7c (cf. Pacholczyk 1970) 
and hence some excess emission (i.e. a local peak) near the cooling frequency Vc may 
be observed. Some of the results presented above find application in modelling the light 
curve and spectrum of GRB 010222 afterglow (Sagar et al 2001, Cowsik et al 2001). 
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